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Abstract. The genus Goniurosaurus (tiger geckos) currently consists of 23 species distributed in China, 
Japan and Vietnam. Several species complexes and recent discoveries of cryptic species pose challenges 
to the species identification, which is crucial to effectively implement the recent listing of the species 
from China and Vietnam in CITES Appendix I and the species from Japan in CITES Appendix III. 
Based on the results of our field work in northern Vietnam and data compiled from literature, we herein 
provide a taxonomic review of the genus Goniurosaurus. Our phylogenetic analyses showed that all 
recorded populations of tiger geckos from Vietnam, which were found to be monophyletic with low 
intra-specific genetic divergences, are assigned to one of the four species: G. catbaensis, G. huuliensis, 
G. lichtenfelderi or G. luii. Both genetic and morphological analyses confirm that the species from China 
and Vietnam can be split into three major groups. Based on the newly collected data, we provide an 
extended morphological description of the Vietnamese species. In addition, we provide an identification 
key for all Goniurosaurus species from China, Japan and Vietnam in order to assist authorities in the 
enforcement of the recent CITES listing. 


Keywords. CITES enforcement, morphology, molecular phylogeny, taxonomy, tiger geckos. 


Ngo H.N., Nguyen H.Q., Tran H.M., Ngo H.T., Le M.D., Gewiss L.R., Schingen-Khan M. van, Nguyen T.Q. & 
Ziegler T. 2021. A morphological and molecular review of the genus Goniurosaurus, including an identification 
key. European Journal of Taxonomy 751: 38-67. https://doi.org/10.5852/ejt.2021.751.1379 


Introduction 


The eublepharid genus Goniurosaurus Grismer, Viets & Boyle, 1999 (tiger geckos) currently consists of 
23 species associated with granitic or karst formations in Southeast and East Asia. This genus shows a 
high level of local endemism and most tiger gecko species have been only recorded from a single locality, 
within unique ecosystems or an isolated archipelago (Orlov et al. 2008; Ziegler et al. 2008; Yang & Chan 
2015; Honda & Ota 2017; Zhou et al. 2018; Ngo et al. 2019a; Qi et al. 2020a, 2020b; Zhu et al. 2020a, 
2020b). The complex topography and variable climatic conditions throughout the distribution range of 
the genus may account for the disjunct distribution and occupation of different ecological niches by 
distinct populations. These factors may have driven natural selection, morphological diversification and 
phylogenetic evolution within this genus (Vitt et al. 1997; Clements et al. 2006; Sexton et al. 2009; Gomes 
et al. 2016; Liang et al. 2018; Qi et al. 2020a, 2020b). According to recent morphological and molecular 
analyses, the genus Goniurosaurus is split into four species groups, namely the G. lichtenfelderi group 
with five species from both mainland and islands in China and Vietnam; the G. kuroiwae group containing 
six species from the Ryukyu Archipelago, Japan; the G. /uii group comprising eight species distributed 
throughout islands and mainland in Vietnam and China; and the G. yingdeensis group consisting of four 
species from China (Nguyen et al. 2009; Nguyen 2011; Wang et a/. 2013; Honda & Ota 2017; Liang 
et al. 2018; Qi et al. 2020a, 2020b; Zhu et al. 2020a, 2020b). However, the systematics of the genus 
Goniurosaurus remains challenging due to ongoing discoveries of further cryptic species and the fact that 
there is not a single genetic marker that covers all 23 species, precluding a complete generic phylogeny of 
the group. Five subspecies of the G. kuroiwae group from Japan were reinstated at full species status and 
anew species, G. sengokui (Honda & Ota, 2017), was recently discovered from the Ryukyu Archipelago 
(Honda & Ota 2017). In China, two taxa, previously identified as G. /uii Grismer, Viets & Boyle, 1999, 
were described as distinct species, namely G. kadoorieorum Yang & Chan, 2015 and G. kwangsiensis 
Yang & Chan, 2015 (Yang & Chan 2015), and a sister species of G. araneus Grismer, Viets & Boyle, 1999, 
namely G. gezhi Zhu, He & Li, 2020, was recently discovered (Zhu et al. 2020a). These three species 
occur in Southwest Guangxi Province (Zhu et al. 2020a). Four additional species, namely G. kwanghua 
Zhu & He, 2020 and G. zhoui Zhou, Wang, Chen & Liang, 2018, belonging to the G. lichtenfelderi group, 
as well as G. gollum Qi, Wang, Grismer, Chen, Lyu & Wang, 2020 and G. varius Qi, Grismer, Lyu, 
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Zhang, Li & Wang, 2020 of the G. yingdeensis group, were recently described from China (Zhou et al. 
2018; Qi et al. 2020a, 2020b; Zhu et al. 2020b). 


To date, five species of Goniurosaurus are reported from Vietnam, viz. G. araneus, G. catbaensis 
Ziegler, Nguyen, Schmitz, Stenke & Résler, 2008, G. huuliensis Orlov, Ryabov, Nguyen, Nguyen & Ho, 
2008, G. lichtenfelderi (Mocquard, 1897) and G. /uii (Grismer et al. 1999; Vu et al. 2006; Orlov et al. 
2008; Ziegler et al. 2008; Nguyen et al. 2009; Nguyen 2011). All species were described based only on 
a few specimens and a small set of diagnostic characters. As such, phenotypic variability among these 
species may lead to the misidentification of taxa. For example, Orlov & Darevsky (1999) described 
a new species (G. murphyi) based on a juvenile specimen, which was subsequently synonymized 
with G. lichtenfelderi by Grismer (2000). Ngo et al. (2016) documented another case with indistinct 
morphological differences between the newly described G. kadoorieorum and G. luii. 


Due to habitat degradation and over-harvesting for the pet trade, wild populations of species of 
Goniurosaurus from China and Vietnam have been subject to severe declines (Stuart et al. 2006; Yang & 
Chan 2015; Ngo et al. 2019b). Therefore, all Chinese and Vietnamese tiger gecko species have recently 
been included in CITES (Convention on International Trade in Endangered Species of Wild Fauna 
and Flora) Appendix II (CITES 2019; Ngo et al. 2019b). To effectively enforce relevant national and 
international regulations and to enable sustainable trade in CITES listed species, a detailed guideline 
for species identification is an essential prerequisite, especially as trade in the morphologically similar 
species from Japan is currently regulated under a lower protection level, namely CITES Appendix II 
(CITES Notification No. 2020/068). A proper identification guidance will help to prevent certain taxa 
from being traded under a wrong name to circumvent trade restrictions. 


In this study, phylogenetic analyses were performed across all recorded populations of Goniurosaurus 
from Vietnam and along with other species from China and Japan using a fragment of the mitochondrial 
Cytochrome b gene. Based on newly collected data from Vietnam, we expanded morphological 
descriptions of each species and evaluated intra- and inter-specific morphological variations among the 
species of Goniurosaurus. We further compared morphological data of Vietnamese tiger geckos with 
those from China and identified characters that can distinguish the Japanese species. We thereby provide 
an identification key for all Goniurosaurus species from China, Japan and Vietnam, in order to assist 
CITES enforcement authorities in species identification. 


Material and methods 
Field surveys 


Field surveys were conducted in June and August 2014, May 2015, July 2016, June and July 2017, 
April 2018, from April to September 2019, as well as in June and July 2020 in Vietnam. Surveys were 
conducted for G. catbaensis in Hai Phong City and Quang Ninh Province, for G. Auuliensis in Lang Son 
Province, for G. lichtenfelderi in Hai Duong, Bac Giang and Quang Ninh provinces, and for G. /uii in Cao 
Bang and Lang Son provinces. Animals were captured by hand and subsequently released at the same 
sites after taking measurements and photos in life. Coordinates of each captured individual were recorded 
using a GPS device (Garmin GPSmap64, WGS84 datum) and can be shared upon request to the authors. 


In addition, several vouchered specimens and small tissue samples of Goniurosaurus deposited in 
collections of the Institute of Ecology and Biological Resources (IEBR), Vietnam National Museum 
of Nature (VNMN) and Hanoi University of Science (HUS), Hanoi, Vietnam, were examined for 
morphological characters and molecular phylogeny. 
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Molecular data and phylogenetic analyses 


All taxa of the genus Goniurosaurus from China, Japan, and Vietnam, for which Cytochrome b data 
were available in GenBank, were included in the study. Three species, Gekko chinensis (Gray, 1842) 
(NC027191), Coleonyx mitratus Peters, 1863 (AB853481) and Holodactylus africanus Boettger, 1893 
(AB853482), were used as outgroups following Honda et al. (2014). 


Small tissue samples were collected from the tail tip of 18 wild specimens from representative localities 
in Vietnam and preserved separately in 70% ethanol (Merck, Germany). The mitochondrial DNA 
Cytochrome b (cytb) gene was selected for DNA sequencing. The total genomic DNA was extracted 
using GeneJet Genomic DNA Purification (ThermoFisher Scientific, Lithuania), following protocols by 
the manufacturer’s instructions. PCR reactions were performed using HotStar Taq Mastermix (Qiagen, 
Germany) to amplify a fragment of approximately 1000 bps of cytb. The primer pair used for this 
study was L14731 (5’°-TGGTCTGAAAAACCATTGTTG-3’; Honda et al. 2014) and GoniRI (5’— 
CTACGGGCTGTCCTCCGATTCAGGTT-:3’; this study). The PCR volume consisted of 21 pl: 2 pl 
of each primer, 5 1 water, 10 pl of Taq mastermix and 1—2 ul DNA template depending on the quality 
of DNA in the final extraction solution. The PCR was performed at 95°C for 15 minutes, followed by 
35 cycles, for 30 s at 95°C, 45 s at 48°C and 60 s at 72°C with a final elongation step for 6 minutes at 
72°C. Negative and positive controls were used for all DNA extractions and PCR reactions. 


PCR products were visualized using electrophoresis through a 1% agarose gel, marker 1kb, 1X TBE, 
stained with ethidium bromide and photographed under UV light. Successful amplifications were purified 
using GeneJet PCR Purification Kit (ThermoFisher Scientific, Lithuania). Cleaned PCR products were 
sent to 1*' Base (Malaysia) for sequencing in both directions. 


The obtained sequences were aligned in Sequencher ver. 5.4 (Gene Codes Corp, Ann Arbor, MI, USA) 
and afterwards aligned using ClustalX ver. 2.1 (Thompson et al. 1997) with default settings. Data were 
analyzed using maximum parsimony (MP) as implemented in PAUP*4.0b10 (Swofford 2001) and 
Bayesian inference (BI) as implemented in MrBayes ver. 3.2 (Ronquist et al. 2012). Settings for these 
analyses followed Le et al. (2006), except that the number of generations in the Bayesian analysis was 
increased to 1 x 10’. For the maximum likelihood (ML) analysis, we used IQ-TREE ver. 1.6.7.1 (Nguyen 
et al. 2015) with a single model and 10000 ultrafast bootstrap replications. The optimal model for 
nucleotide evolution was set to TPM1uf+I+G as selected by Jmodeltest ver. 2.1.4 (Posada et al. 1998). 
As the phylogenetic relationships between species and species groups were well resolved and virtually 
every important node received high statistical support from all analyses, we opted not to partition our 
data by codon positions. 


For the Bayesian analysis, we used the optimal model determined by Jmodeltest with parameters 
estimated by MrBayes ver. 3.2.7. Four Markov chains, one cold and three heated utilizing default 
heating values, were sampled every 1000 generations. Log-likelihood scores of sample points were 
plotted against generation time to detect stationarity of the Markov chains. The burn-in value was set 
to 50 in the BI analysis, as -InL scores reached stationarity after 50 000 generations in both runs. Nodal 
support was evaluated using Bootstrap replication (BP) as estimated in PAUP*4.0b10 and ultrafast BP in 
IQ-TREE ver. 1.6.7.1 and posterior probability (PP) in MrBayes ver. 3.2. BP > 70 (Hillis & Bull 1993) 
and ultrafast BP (UBP) and PP > 95% are regarded as strong support for a clade (Ronquist et al. 2012; 
Nguyen et al. 2015). The uncorrected pairwise distance (p) were calculated in PAUP*4.0b10. 


Morphological analyses 


A total of 486 live individuals and 54 museum specimens of four species from Vietnam were examined 
for morphological data, comprising 194 individuals of G. catbaensis (21 juveniles, 93 females and 
80 males), 80 individuals of G. huuliensis (02 juveniles, 46 females and 32 males), 178 individuals of 
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G. lichtenfelderi (14 juveniles, 72 females and 92 males), and 88 individuals of G. /uii (11 juveniles, 
43 females and 34 males). Detailed descriptions of each species were based on newly collected data 
combined with previous literature for the morphological variation, such as for G. catbaensis from 
Ziegler et al. (2008) and Nguyen (2011), for G. huuliensis from Orlov et al. (2008) and Nguyen (2011), 
for G. lichtenfelderi from Grismer (2000), Grismer et a/. (2002) and Nguyen (2011), and for G. /uii from 
Grismer et al. (1999), Vu et al. (2006) and Nguyen (2011). No specimen of G. araneus was investigated 
in this study, however, we included morphological data of eight specimens from the type locality in Cao 
Bang Province, Vietnam and Guangxi Province, China from the literature (Grismer et al. 1999; Chen 
et al. 2014). The sex of each collected specimen was determined based on the presence (in males) or 
absence (in females) of large swollen hemipenial bulges. Lizards were categorized into two age classes 
based on the snout-vent length (G. lichtenfelderi: SVL < 80 mm=juveniles, SVL > 80 mm = adults, 
while three other species of Goniurosaurus were sorted with SVL < 85 mm=juveniles, SVL > 85 mm = 
adults). 


Measurements were taken with dial calipers to the nearest 0.1 mm at the right side of each individual. 


Abbreviations 

AD = _ diameter of auditory meatus 

AG = axilla to groin length, from posterior edge of forelimb insertion to anterior edge of hindlimb 
insertion 

BH = maximum body height, from top of dorsal body to belly 

BW = maximum body width, greatest width of torso, taken at level of midbody 

CH = cheek height, from posterior edge of labial to top of head at parietal region 

ED = diameter of eye, greatest diameter of orbit 

EE = eye to ear distance, from posterior margin of eye to posterior margin of ear 

FLL = forelimb length, from axilla to the tip of the fourth finger 

HH = maximum head height 

HL = _ head length, from the tip of snout to posterior edge of occiput 


HLL = hindlimb length, from groin to the tip of the fourth toe 
HW = maximum head width 


IO1 = interorbital distance, distance between anteriormost points of eyes 

102 = _ interorbital distance, distance between posteriormost points of eyes 

ML = mouth length, from tip of snout to last posterior labial edge 

MW = _ mouth width, distance between last posterior labial edges on each side 

ND =  supranasal distance, distance between nares 

SE = snout to eye distance, measured from tip of snout to anteriormost point of eye 
SVL = snout-vent length, from tip of snout to vent 

TaL = tail length, from vent to tip of tail 

WT = maximum tail width 


Scale counts 


CIL = eyelid fringe scales or ciliaria 

DTR = dorsal tubercle rows at midbody 

GP = gular scales bordering the postmentals 

GST = granular scales surrounding dorsal tubercles 
IFL = infralabials 

IN = postrostrals or internasals 


LD1 = subdigital lamellae under the first finger 
LD4 = subdigital lamellae under the fourth finger 
LT1 =  subdigital lamellae under the first toe 
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LT4 = subdigital lamellae under the fourth toe 

MB_ = scales around midbody 

N = nasal scales surrounding nare 

P-IN = post-internasals 

PAT = postcloacal tubercles 

PM = postmentals 

PO = preorbital scales 

PP = precloacal pores 

SPL = _ supralabials 

TL = paravertebral tubercles between limb insertions 


Statistical analyses were performed by using the software environment R.3.1.2 (RStudio Team 2018). 
Shapiro- Wilk’s test was used to test the assumption of normality. Kruskal-Wallis’s test was performed to 
determine differences concerning the SVL of sex and age classes among species, as well as the extracted 
PC1 and PC2 values (mentioned below) of a principal component analysis among four tiger gecko 
species in Vietnam except of G. araneus. For both tests we applied p=0.05. The SVL variable was 
excluded from subsequent analyses due to high collinearity with other dimensions. We further performed 
a Principal Component Analysis (PCA) of log,,-transformed raw data of 18 remaining morphometric 
characters (except TaL and WT variables, due to many regenerated tails), using the packages “factoextra” 
and “FactoMinerR” to detect variances among the four selected Vietnamese species via contribution 
percentage of the PC1 and other PCs scores (Kassambara 2017). Morphometric variation of each species 
was illustrated by representatively clustered ellipse spaces with different coded-color within a spatial 
coordinate of the first two most important dimensions (PC1 and PC2) in the PCA analysis, to visually 
evaluate the overlap among the four species. In addition, we identified the overall difference in meristic 
characters among the 17 tiger gecko species native to China and Vietnam assigned to three groups 
(G. lichtenfelderi, G. luii and G. yingdeensis) by using a Multiple Correspondence Analysis (MCA). 
The meristic variation among the three groups was also visualized by convex ellipses within a spatial 
coordinate of the first two most important dimensions (Dim! and Dim2) in the MCA analysis. 


Results 


Molecular phylogeny 


The final matrix consisted of 949 bp aligned characters with 48 sequences of 14 ingroup and three outgroup 
taxa. The alignment contained no gaps. In total, 404 characters were found to be parsimony informative. 
The MP analysis produced 624 most parsimonious trees with 1404 steps (CI=0.54, RI = 0.85). Our 
phylogenetic analyses recovered a generally similar topology to those reported by Liang et al. (2018) 
and Zhu et al. (2020a, 2020b). Specifically, the Goniurosaurus lichtenfelderi group and the G. luii group 
form a monophyletic group with high nodal values from all analyses and the G. kuroiwae group is sister 
to the other two. In addition, each species group was also strongly corroborated as monophyletic from all 
analyses (Fig. 1). Within each species group, the results supported by this study are more similar to those 
generated by Liang et al. (2018), because Zhu et al. (2020a) used a different mitochondrial marker, the 
16S ribosomal RNA gene. Nonetheless, we found several noted discrepancies between the two studies. 
Specifically, our phylogenetic analyses showed that G. kuroiwae (Namiye, 1912) is polyphyletic with high 
nodal support for each of the clades. In addition, G. /uii was strongly corroborated as monophyletic in our 
Bayesian and MP analyses (BP=93%; PP=98%), but only weakly recovered as monophyletic in Liang 
et al. (2018) (BP=67%; PP=64%). Goniurosaurus lichtenfelderi was recovered as monophyletic with 
perfect nodal support from all analyses in our study, but G. hainanensis Barbour, 1908 is polyphyletic. 
On the other hand, both species are not monophyletic in Liang et al. (2018). The positions of G. zhoui 
and G. bawanglingensis Grismer, Haitao, Orlov & Anajeva, 2002 are interchanged in the trees supported 
by the two studies with strong nodal support from Liang et a/. (2018) and weak corroboration from 
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ours. All surveyed sub-populations of Goniurosaurus from Vietnam were placed in four separate clades 
corresponding to four distinct species, G. catbaensis, G. huuliensis, G. luii and G. lichtenfelderi with 
significantly statistical support values, except for the ML analysis of G. /uii (UBP=88%). Intra-specific 
genetic divergences of these species were relatively low: G. catbaensis (0.7—1.6%), G. huuliensis (0.4— 
1.2%), G. lichtenfelderi (0O-0.7%), and G. luii (0—1.9%) (Fig. 1; Supp. file 1: Table S1). 


Taxonomic accounts 


Class Squamata Oppel, 1811 
Order Gekkota Cuvier, 1817 
Family Eublepharidae Boulenger, 1883 
Genus Goniurosaurus Grismer, Viets & Boyle, 1999 


Goniurosaurus araneus Grismer, Viets & Boyle, 1999 


Diagnosis 


Body splayed and gracile, SVL 108.9-124.0 mm; external nares bordered by 6—8 nasal scales; supraorbital 
region with a row of enlarged tubercles; outer surface of upper eyelid composed of small granular 
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Fig. 1. Phylogram containing Goniurosaurus Grismer, Viets & Boyle, 1999 representatives from China, 
Japan and Vietnam based on the Bayesian analysis. Numbers above and below branches are MP/ML 
ultrafast bootstrap values and combined Bayesian posterior probabilities, respectively. Dash represents 
<50% value. 


44 


NGO HN. et al., Review of the genus Goniurosaurus 


scales, one-half the size of those on top of head; internasal single; supralabials 8—10; infralabials 8 or 9; 
preorbital scales 13-18; eyelid fringe scales 52-67; postmentals 4-6; dorsal body scales elongate; 
paravertebral tubercles 29-38; scales around midbody 129-147; scales surrounding dorsal tubercles 10— 
14; axillary pockets deep; subdigital lamellae under fourth toe 23—25; precloacal pores in males 18—23; 
iris dark brown; dorsal ground color of head, body and limbs immaculate dull yellow-gray; nuchal loop 
posteriorly protracted (in V-shape); dorsal body bands between limb insertions 3, wide, edged anteriorly 
and posteriorly by wide dark brown bands; light band on tail base in width of 15 or 16 granular scales; 
ground color of tail black, caudal bands 5, white, completed ventrally; ventral surface of head, body and 
limbs dull white and immaculate (Grismer et a/. 1999; Nguyen 2011; Chen et al. 2014). 


Remarks 


During our field work in northern Vietnam, including the type locality in Cao Bang Province, no specimen 
of Goniurosaurus araneus was seen. Therefore, the diagnosis was solely based on the descriptions of 
Grismer et al. (1999) and Chen et al. (2014). 


Goniurosaurus catbaensis Ziegler, Nguyen, Schmitz, Stenke & Résler, 2008 
Fig. 2 


Diagnosis 

Body splayed and gracile, SVL 89.3—125.3 mm; external nares bordered by 6-8 nasal scales; supraorbital 
region with a row of enlarged tubercles; outer surface of upper eyelid composed of granular scales, about 
the same size of those on top of head and with a row of 6-10 enlarged tubercles; internasals absent; 
supralabials 8-11; infralabials 7-10; eyelid fringe scales 45-56; postmentals 2—5; gular region below 
lower jaws without enlarged tubercles; paravertebral tubercles 31-38; scale rows around midbody 112— 
127, granular scales surrounding tubercles 8—11; axillary pockets deep; subdigital lamellae under fourth 
toe 22—25; precloacal pores in males 16—23; iris orange-brown; dorsal ground color of head, body and 
limbs gray-brown to pale brown and mottled with dark brown blotches; nuchal loop thin, posteriorly 
protracted (in V-shape); dorsal body bands between limb insertions 3-4, thin, yellow, without dark 
spots; light band on tail base in width of 8—9 scales; ground color of tail black, caudal bands 5, white, 
completed ventrally; ventral surface of head, body and limbs dull white and immaculate, gular region 
with brown spots (modified after Ziegler et al. 2008; Nguyen 2011). 


Description (Supp. file 1: Table S2) 


Body splayed and gracile, adult males: SVL92.9-125.3 mm (mean+ SE: 112.3+0.8 mm, n= 80), TaL 7.2— 
97.9 mm (75.2+1.9 mm); adult females: SVL 89.3-122.1 mm (111.8+0.8 mm, n=93), TaL 17.5— 
98.3 mm (70.3+1.7 mm); juveniles: SVL 53.4-78.8 mm (68.8+1.8 mm, n=21), TaL 11.8-61.2 mm 
(48.1+2.6 mm) (Supp. file 1: Table S2); head triangular, wider than neck, covered by uniform granular 
scales interspersed with tubercles in temporal and occipital regions; scales on rostrum slightly larger 
and flatter; enlarged supraorbital tubercles in a conspicuous row; middorsal portion of rostral partially 
sutured dorsomedially, bordered laterally by first supralabial on each side, dorsolaterally by prenasal on 
each side, and dorsally by two supranasals; internasal (postrostral) scales absent; external nares bordered 
by 6-8 nasals: anteriorly by prenasal and supranasal, dorsally by supranasal, posteriorly by two slightly 
enlarged postnasals and | or 2 smaller granular scales, and ventrally by prenasal; prenasals with long 
recurved ventral portion; supranasals triangular, meeting in midline behind rostral suture; preorbital 
scales 10-13; supralabials 8-11; infralabials 7-10; eyes relatively large, pupils vertical; eyelid fringe 
scales 45—56, those of upper eyelid slightly enlarged; outer surface of upper eyelid composed of granular 
scales of about the same size of those on top of head, including a row of 6—10 enlarged tubercles; fold of 
skin originating from suborbital region extends posteroventrally across angle of jaw; external auditory 
meatus elliptical with long axis directed dorsoventrally; tympanum deeply recessed; mental triangular, 
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Fig. 2. Goniurosaurus catbaensis Ziegler, Nguyen, Schmitz, Stenke & Résler, 2008. A. Juvenile 
(SVL=55 mm). B. Juvenile (SVL=69 mm). C. Juvenile (SVL=75 mm). D. Adult (SVL=91 mm). 
E. Adult (SVL=99 mm). F. Adult (SVL=105 mm). G. Adult (SVL=116 mm). H. Adult (SVL=125 mm). 
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bordered laterally by first infralabial on each side and posteriorly by 2—5 postmentals; postmentals 
bordered by 6-10 gular scales; gular region below lower jaws without enlarged tubercles; gular scales 
juxtaposed and granular, abruptly grading posteriorly into flat hexagonal pectoral scales and even larger 
ventral scales. 


Neck narrower than body, covered with uniform granular scales interspersed with several sharply pointed 
conical tubercles on nape; tubercles on flanks conical, those of vertebral region somewhat lower in profile; 
dorsal body tubercles surrounded by 9-11 granular scales; dorsal tubercle rows at midbody 19-24; 
paravertebral tubercles between limb insertions 3 1—38, distinct vertebral row of tubercles absent; scales 
around midbody 112-127; larger ventral scales grade abruptly into smaller granular scales immediately 
anterior to vent at level of preanal pores; males with 18—23 precloacal pores in a transverse continuous 
series; region posterior to vent covered by flat juxtaposed scales and great hemipenial bulges, containing 
2 or 3 enlarged postcloacal tubercles laterally on each side at level of vent; tail long and thin, thick at 
base, anteriorly with whorls; dorsal tail scales flat, smooth, up to 1.5 times the size of dorsal scales, 
arranged in more or less regular transverse rows; light band on tail base in width of 8—9 granular scales 
and with 7-8 tubercles in a transversal series; subcaudals larger than those on dorsal surface of tail. 


Limbs relatively long and thin, covered dorsally with granular scales interspersed with several tubercles 
and ventrally with flat juxtaposed to subimbricate scales; dorsal granular scales grade into slightly 
flattened subimbricate scales on top of pes and manus; hind limbs slightly larger than forelimbs; larger 
granular scales on ventral surface of pes and manus; axillary pockets deep; subdigital lamellae wide, 
9-11 under first finger, 18—21 under fourth finger, 9-12 under first toe, 22—25 under fourth toe; digits 
laterally compressed, increasing in length from first to fourth, fifth shorter than fourth (modified after 
Ziegler et al. 2008; Nguyen 2011). 


Coloration in life (Fig. 2) 


Dorsal ground color of head, body and limbs grey brown (in animals with SVL 53.4-105 mm) and 
blotches chestnut brown (in animals with SVL 110—125.3 mm), juveniles without small blotches, adults 
mottled with few circular blotches on body and limbs and long dark brown blotches on head; few yellow 
conical tubercles on franks of neck and body, and limbs in adults; iris light orange or red brown; five 
bands on the dorsal ground, thin, immaculate without dark spots, cream in juveniles and light orange 
or yellow in adults, all edged anteriorly and posteriorly by thin dark brown bands, including one thin 
nuchal loop extending from posterior corners of eyes and posteriorly protracted (in V-shape), three body 
bands between limb insertions, and another one on tail base; ground color of tail dark brown, and grey 
brown at mottled tail base; 3-5 immaculate white caudal bands, edged anteriorly and posteriorly in 
black; ventral surface of head, body and limbs dull white, juveniles absolutely immaculate, but adults 
with few dark spots on limbs, weak brown lateral spotting in gular region, venter and limbs. 


Goniurosaurus huuliensis Orlov, Ryabov, Nguyen, Nguyen & Ho, 2008 
Fig. 3 


Diagnosis 

Body splayed and gracile, SVL 97.2—134.6 mm; external nares bordered by 6-8 nasal scales; supraorbital 
region with a row of enlarged tubercles; outer surface of upper eyelid composed of granular scales, about 
one half the size of those on top of head and without enlarged tubercles; internasal | or 2 (rarely absent); 
supralabials 9-12; infralabials 9-12; preorbital scales 14-20; eyelid fringe scales 51-59; postmentals 
2-4; gular region below lower jaws with enlarged tubercles; paravertebral tubercles 31—37; scale rows 
around midbody 118-130, granular scales surrounding tubercles 11—13; axillary pockets deep; subdigital 
lamellae under fourth toe 21—25; precloacal pores in males 25—30; iris reddish brown; dorsal ground 
color of head, body and limbs dark brown, without small dark brown blotches (dark blotches present 
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only on lower zone of flanks); nuchal loop thin, posteriorly protracted (in V-shape); dorsal body bands 
between limb insertions 3, thin, immaculate yellow; gular region with brown spots (modified after Orlov 
et al. 2008; Nguyen 2011). 


Description (Supp. file 1: Table S2) 


Body splayed and gracile, adult males SVL 97.2-132.2 mm (mean+SE: 118.9+1.4 mm, n=32), 
TaL 36.7-108.6 mm (72.7+3.7 mm), adult females SVL 97.4-134.6 mm (121.1+1.2 mm, n=46), 


Fig. 3. Goniurosaurus huuliensis Orlov, Ryabov, Nguyen, Nguyen & Ho, 2008. A. Juvenile 
(SVL=74.4 mm). B. Adult (SVL=97 mm). C. Adult (SVL=108 mm). D. Adult (SVL=116 mm). 
E. Adult (SVL= 125 mm). F. Adult (SVL=132 mm). 
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TaL 48.4-106.7 mm (78.2+2.7 mm), juveniles SVL 73.5-74.4 mm (n=2), TaL 62.5-62.7 mm 
(Supp. file 1: Table S2); head triangular, wider than neck, covered by uniform granular scales interspersed 
with tubercles in temporal and occipital regions; scales on rostrum slightly larger and flatter; enlarged 
supraorbital tubercles in a conspicuous row; middorsal portion of rostral partially sutured dorsomedially, 
bordered laterally by first supralabial on each side, dorsolaterally by prenasal on each side, and dorsally 
by internasal and two supranasals; internasal 1-3 (rarely absent); external nares bordered by 5—7 nasals 
(prenasal, anterior and posterior supranasals, 2 slightly enlarged postnasal, and | or 2 granular scales); 
preorbital scales 14-20; supralabials 9-12, grading into granular scales posteriorly; infralabials 9-12; 
eyes large, pupils vertical; eyelid fringe scales 51-59, those of upper eyelid slightly enlarged; outer 
surface of upper eyelid composed of granular scales of about one half the size of those on top of head, 
without enlarged tubercles; fold of skin originating from suborbital region extends posteroventrally 
across angle of jaw; external auditory meatus elliptical with long axis directed dorsoventrally; tympanum 
deeply recessed; mental triangular, bordered laterally by first infralabial on each side and posteriorly 
by 2-4 postmentals; postmentals bordered by 7—10 gular scales; gular region below lower jaws with 
enlarged tubercles; gular scales juxtaposed and granular, grading posteriorly into flat hexagonal pectoral 
scales and even larger ventral scales. 


Neck narrower than body, covered with uniform granular scales interspersed with several sharply 
pointed conical tubercles on nape; tubercles on flanks conical, those of vertebral region somewhat 
flatter; dorsal body tubercles surrounded by 11-13 granular scales; dorsal tubercle rows at midbody 
19-24; paravertebral tubercles between limb insertions 3 1-37, distinct vertebral row of tubercles absent; 
scales around midbody 118—130; ventral scales large; males with 25—30 precloacal pores in continuous 
series, females without precloacal pores (but pitted scales present); region posterior to vent covered by 
flat juxtaposed scales and greatly swollen, containing | or 2 enlarged tubercles on each side at level of 
vent; tail thick at base; light band on tail base in width of 9-12 granular scales and with 8—10 tubercles 
in transversal series. 


Limbs relatively long and thin, covered dorsally with granular scales interspersed with several tubercles 
and ventrally with flat juxtaposed to subimbricate scales; dorsal granular scales grade into slightly 
flattened subimbricate scales on top of pes and manus; hind limbs larger than forelimbs; axillary pockets 
deep; subdigital lamellae wide, 10 or 11 under first finger, 18—21 under fourth finger, 11 or 12 under first 
toe, 21-25 under fourth toe; digits laterally compressed, increasing in length from first to fourth, fifth 
shorter than fourth (modified after Orlov et al. 2008; Nguyen 2011). 


Coloration in life (Fig. 3) 


Dorsal ground color of head, body and limbs signal brown (juveniles), grey-brown (young adults, SVL: 
97.2-116 mm) and signal brown (adults), without small dark brown blotches (dark blotches present only 
on lower zone of flanks); dull white tubercles on dorsal body, limbs in juveniles, few orange tubercles on 
limbs in adults; iris orange or red brown; five bands on the dorsal ground, thin, immaculate without dark 
spots, slight bisque in juveniles, orange brown or yellow in adults, all edged anteriorly and posteriorly 
by thin dark brown bands, including one thin nuchal loop extending from posterior corners of eyes and 
posteriorly protracted (in V-shape), three body bands between limb insertions, and another one on tail 
base; ground color of tail dark brown, and signal brown at mottled tail base; 3-6 immaculate white 
caudal bands; ventral surfaces of head, body and limbs dull white and immaculate except for few dark 
brown spots on margin regions of gular and limbs. 
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Goniurosaurus luii Grismer, Viets & Boyle, 1999 
Fig. 4 


Diagnosis 

Body splayed and gracile, SVL 86.5—126.5 mm; external nares bordered by 5—8 nasal scales; supraorbital 
region with a row of enlarged tubercles; outer surface of upper eyelid composed of granular scales, 
about one half the size of those on top of head and without enlarged tubercles; internasals 1 or 2; 
supralabials 8—12; infralabials 8-11; preorbital scales 13-16; eyelid fringe scales 46-61; postmentals 
2-6; gular region below lower jaws with enlarged tubercles; paravertebral tubercles 29-38; scale rows 
around midbody 119-— 144, granular scales surrounding tubercles 9—14; axillary pockets deep; subdigital 
lamellae under fourth toe 20—26; precloacal pores in males 23-32; iris brown or bright orange; dorsal 
ground color of head, body and limbs pale brown to grey brown, mottled with small dark brown blotched; 
nuchal loop thin, posteriorly protracted (in V-shape); dorsal body bands between limb insertions 3, thin, 
immaculate yellow; gular region, belly, and ventral surface of limbs with brown spots (Grismer et al. 
1999; Vu et al. 2006; Nguyen 2011). 


Description (Supp. file 1: Table S2) 


Body splayed and gracile, males: SVL 88.8-123.0 mm (mean+SE: 109.34+1.4 mm, n=34), 
TaL 2.8-96.8 mm (72.5+3.8 mm); adult females: SVL 86.5—126.5 mm (112.4+1.3 mm, n=43), TaL 
43.3-102.2 mm (72.2+2.1 mm); juveniles: SVL 55.0-84.8 mm (74.8+3.1 mm, n=11), TaL 44.2— 
68.9 (58.742.7 mm) (Supp. file 1: Table S2); head triangular, wider than neck, covered by uniform 
granular scales interspersed with tubercles in temporal and occipital regions; scales on rostrum slightly 
larger and flatter; enlarged supraorbital tubercles in a conspicuous row; middorsal portion of rostral 
partially sutured dorsomedially, bordered laterally by first supralabial on each side, dorsolaterally by 
prenasal on each side, and dorsally by 1 or 2 internasal and two supranasals; internasals | or 2; external 
nares bordered by 5-8 nasal scales (prenasal, anterior and posterior supranasals, 2 slightly enlarged 
postnasal, and 1-3 granular scales); preorbital scales 13-16; supralabials 8-12, grading into granular 
scales posteriorly; infralabials 8-11; eyes large, pupils vertical; eyelid fringe scales 46-56, those of 
upper eyelid slightly enlarged; outer surface of upper eyelid composed of granular scales, about one half 
the size of those on top of head and without enlarged tubercles; fold of skin originating in suborbital 
region extends posteroventrally across angle of jaw; external auditory meatus elliptical with long 
axis directed dorsoventrally; tympanum deeply recessed; mental triangular, bordered laterally by first 
infralabial on each side and posteriorly by 2—6 postmentals; postmentals bordered by 6—11 gular scales; 
gular region below lower jaws with enlarged tubercles; gular scales juxtaposed and granular, abruptly 
grading posteriorly into flat hexagonal pectoral scales and even larger ventral scales. 


Neck narrower than body, covered with uniform granular scales interspersed with several sharply 
pointed conical tubercles on nape; tubercles on flanks conical, those of vertebral region somewhat more 
flat; dorsal body tubercles surrounded by 11—13 granular scales; dorsal tubercle rows at midbody 20—24; 
paravertebral tubercles between limb insertions 29-38, distinct vertebral row of tubercles absent; scale 
rows around midbody 119-144; ventral scales large; males with 24-32 precloacal pores in a transverse 
continuous series, females without precloacal pores (but pitted scales present); region posterior to vent 
covered by flat juxtaposed scales and greatly swollen, containing 1—3 enlarged tubercles on each side 
at level of vent; tail thick at base, light band on tail base in width of 9-13 granular scales with 9-10 
tubercles in transversal series; ventral caudals of tail base larger than dorsal caudals. 


Limbs relatively long and thin, covered dorsally with granular scales interspersed with several tubercles 
and ventrally with flat juxtaposed to subimbricate scales; dorsal granular scales grade into slightly 
flattened subimbricate scales on top of pes and manus; hind limbs larger than forelimbs; axillary pockets 
deep; subdigital lamellae wide, 9-12 under first finger, 17—22 under fourth finger, 10—12 under first toe, 
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Fig. 4. Goniurosaurus luii Grismer, Viets & Boyle, 1999. A. Juvenile (SVL=61 mm). B. Juvenile 
(SVL=70 mm). C. Juvenile (SVL=80 mm). D. Adult (SVL=89 mm). E. Adult (SVL=96 mm). F. Adult 
(SVL70= 104 mm). G. Adult (SVL=116 mm). H. Adult (SVL=126 mm). 
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20-26 under fourth toe; digits laterally compressed, increasing in length from first to fourth, fifth shorter 
than fourth (Grismer et al. 1999; Vu et al. 2006; Nguyen 2011). 


Coloration in life (Fig. 4) 


Dorsal ground color of head, body and limbs pale brown to grey brown, mottled with small dark 
brown circle blotches, juveniles without small dark brown blotches (except juveniles SVL > 72 mm 
initially present on head); iris red brown or bright orange; five bands on dorsal ground, thin, immaculate 
without dark spots, slight bisque in juveniles, orange brown or yellow in adults, all edged anteriorly 
and posteriorly by thin dark brown bands, including one thin nuchal loop extending from posterior 
corners of eyes and posteriorly protracted (in V-shape), three body bands between limb insertions, and 
another one on tail base; ground color of tail nearly solid black, and grey brown at mottled tail base; 3-6 
immaculate white caudal bands; ventral surfaces of head, body and limbs dull white with a few dark 
brown spots in gular region, on belly and limbs. 


Goniurosaurus lichtenfelderi (Mocquard, 1897) 
Fig. 5 


Diagnosis 

Body robust; SVL 80.6-113.5 mm; external nares bordered by 5—10 nasal scales; supraorbital region 
with a row of slightly enlarged tubercles; outer surface of upper eyelid composed of granular scales, 
about one-half the size of those on top of head and without enlarged tubercles; internasals 1—5 (rarely 
1:2; 2:2 or 2:3); supralabials 7-10; infralabials 6—9; preorbital scales 12—19; eyelid fringe scales 43-58; 
postmentals 2—6; paravertebral tubercles 22—33; scale rows around midbody 117-131, granular scales 
surrounding tubercles 10-13; axillary pockets shallow; subdigital lamellae under fourth toe 17-21; 
precloacal pores in males 25—33, in females 17—21; dorsal ground color of head, body and limbs dark 
purple-brown, without small dark brown blotches; transverse body bands 4, nuchal loop thin, posteriorly 
rounded, in U-shape; dorsal body bands between limb insertions 2, thin, light yellow; gular region and 
without dark spots; ventral surfaces of head, body and limbs dull white with a few dark dots on margin 
regions (modified after Grismer 2000; Grismer et al. 2002; Nguyen 2011). 


Description (Supp. file 1: Table S2) 


Body robust, adult males: SVL 80.6—-113.5 mm (mean+SE: 97.840.7 mm, n=92), TaL 6.2— 
84.8 mm (61.5+1.7 mm); adult females: SVL 81.0-105.5 mm (96.040.7 mm, n=72), TaL 22.7— 
81.2 mm (58.4+1.4 mm); juveniles: SVL 41.1—77.1 mm (65.442.9 mm, n=14), TaL 20.9-64.7 mm 
(48.3+3.4 mm) (Supp. file 1: Table S2); head triangular, wider than neck, covered by uniform granular 
scales interspersed with tubercles in temporal and occipital regions; scales on rostrum slightly larger 
and flatter; enlarged supraorbital tubercles in a conspicuous row; middorsal portion of rostral partially 
sutured dorsomedially, bordered laterally by first supralabial on each side, dorsolaterally by prenasal 
on each side, and dorsally by 1 or 2 internasals and two supranasals; internasals 1—5 (rarely 1:2; 2:2 
or 2:3); external nares bordered by 5—10 nasals; preorbital scales 12-18; supralabials 7-10, grading 
into granular scales posteriorly; infralabials 6—9; eyes large, pupils vertical; eyelid fringe scales 47-58, 
those of upper eyelid slightly enlarged; outer surface of upper eyelid composed of granular scales of 
about one-half the size of those on top of head, without enlarged tubercles; fold of skin originating in 
the suborbital region extends posteroventrally across angle of jaw; external auditory meatus elliptical; 
tympanum deeply recessed; mental triangular, bordered laterally by first infralabial on each side and 
posteriorly by 2—5 postmentals; postmentals bordered by 7-10 gular scales; gular region below lower 
jaws without enlarged tubercles; gular scales juxtaposed and granular, abruptly grading posteriorly into 
flat hexagonal scales and even larger ventral scales. 
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Fig. 5. Goniurosaurus lichtenfelderi (Mocquard, 1897). A. Juvenile (SVL=55 mm). B. Juvenile 
(SVL=63 mm). C. Juvenile (SVL=72 mm). D. Adult (SVL=80 mm). E. Adult (SVL=91 mm). F. Adult 
(SVL=97 mm). G. Adult (SVL=105 mm). H. Adult (SVL=113 mm). 
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Neck narrower than body, covered with uniform granular scales interspersed with several sharply conical 
tubercles on nape; tubercles on flanks conical, those of vertebral region somewhat more flat; dorsal body 
tubercles surrounded by 10-13 granular scales; dorsal tubercle rows at midbody 15-22; paravertebral 
tubercles between limb insertions 22-33, distinct vertebral row of tubercles absent; scales around 
midbody 117-130; ventral scales large; males with 24—33 precloacal pores in a transverse continuous 
series, females without distinct precloacal pores (but 25—33 pitted scales present); region posterior to 
vent covered by flat juxtaposed scales and greatly swollen, containing | (rarely 2) enlarged tubercles on 
each side at the level of vent; tail thick at base; light band on tail base in width of 4—7 scales and with 
6-10 tubercles in a transversal series. 


Limbs robust, covered dorsally with granular scales interspersed with several tubercles and ventrally with 
flat juxtaposed to subimbricate scales; dorsal granular scales grade into slightly flattened subimbricate 
scales on top of pes and manus; hind limbs larger than forelimbs; axillary pockets shallow; subdigital 
lamellae wide, 6—9 under first finger, 13-17 under fourth finger, 6-10 under first toe, 16—20 under 
fourth toe; digits laterally compressed, increasing in length from first to fourth, fifth shorter than fourth 
(modified after Grismer 2000; Grismer et al. 2002; Nguyen 2011). 


Coloration in life (Fig. 5) 


Dorsal ground color of head, body and limbs grey and grey brown in juveniles without blotches and 
chestnut brown in adults with dark brown blotches (in animals with SVL > 80 mm); iris orange or 
redish brown; four bands on the dorsal ground, thin, yellow in both juveniles and adults, immaculate 
without dark spots (few dark brown spots present in some animals), all edged anteriorly and posteriorly 
by thin dark brown bands, including one thin nuchal loop extending from posterior corners of eyes and 
posteriorly rounded (in U-shape), two body bands between limb insertions, and another one on tail base; 
ground color of tail dark brown, and grey brown at mottled tail base; 3—5 immaculate white caudal bands 
except first band slightly yellow, and some animals with regenerated tail present immaculate grey brown 
without white caudal bands; ventral surfaces of head, body and limbs dull white with a few dark dots 
on margin regions of belly, limbs, gular region, but immaculate dull white in juveniles (SVL < 65 mm) 
without dark spots. 


Morphological comparisons 


Our morphological comparisons of the tiger gecko species (except for G. araneus) from Vietnam 
revealed an overall significant difference in the snout-vent length (SVL). As the result, G. huuliensis had 
the largest SVL, while G. lichtenfelderi had the shortest SVL in both adult males and females, compared 
to the other tiger geckos (P < 0.05). However, the SVL of G. catbaensis was not significantly different 
from that of G. /uii (P > 0.05, Supp. file 2: Fig. $1). 


The PCA of 18 selected morphometric characters estimated the first (PC1) and second principal (PC2) 
components to explain 82.5% (74.5% and 8.0%, respectively) of the variance among the four investigated 
species of Goniurosaurus in Vietnam (Supp. file 2: Fig. S2). The extracted PC1 and PC2 scores of the 
PCA differ significantly among the four species (Kruskal-Wallis’s test, P < 0.05), of which the values 
and morphometric spaces of G. lichtenfelderi are not concordant with the three remainders (Fig. S2). 
Three morphological factors, namely head length (HL), head width (HW), and mouth width (MW) 
highly account for the overall difference in the PC1 score, while the PC2 score is mainly explained by 
the body height (BH) and body width (BW) (Supp. file 2: Fig. $2). 


Regarding the comparisons of the 17 species of Goniurosaurus from China and Vietnam, the multiple 
correspondence analysis (MCA) clustered them separately into three groups (G. lichtenfelderi group, 
G. luii group and G. yingdeensis group) represented with different meristic spaces on the spatial 
coordinate of Dim 1 and Dim 2 (Fig. 6A, Table 1). The characters of fourth toe (LT4), body limb shape 
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Fig. 6. A. Multiple correspondence analysis (MCA) on meristic variation among 17 recorded species of 
Goniurosaurus Grismer, Viets & Boyle, 1999 in China and Vietnam. B. Correlation between meristic 
variables and principal dimensions (Dim1 and Dim2). 


and body bands highly account for the Dim1 score, whereas the characters of precloacal pores (PP) 
and granular scales surrounding dorsal tubercles (GST) are the most correlated with the Dim 2 score 
(Fig. 6B). 


Key to the species of Goniurosaurus (Fig. 7, Table 1) 


Modified from Grismer et al. (1999, 2002), Vu et al. (2006), Orlov et al. (2008), Ziegler et al. (2008), 
Wang et al. (2010, 2013, 2014), Nguyen (2011), Chen et al. (2014), Yang & Chan (2015), Honda & Ota 
(2017), Zhou et al. (2018, 2020a, 2020b), Qi et al. (2020a, 2020b) and the present study. 


1. 


Precloacal pores in males present, claws are sheathed by scales .........cceseseseesceseceeeeeeeeeeeeneeneeenees 7 
Precloacal pores absent, unsheathed claws (G. kuroiwde QTOUP) ......:escceseesseseeeeceeceetceeeeeseeaeeneeenees 2 
Yellow brown to gold iris; a single scale at the base of each digit which is occasionally slightly 


CM ALO visadsceaduiadesdeesseatiles aves ishavead axsdediseeasseayedes Gabavee cash aundeeauiaiaants G. yamashinae (Okada, 1936) 
Blood-red iris; one to three enlarged scales at the base of each digit on the limbs ............ eee 3 
Adults without mterspace Motte si cecesisasecassWiacesteessestveseves veaweabehernnthndeatsseictsinsscnscatenvlnvan dunes. 4 
Adults with interspace: MOttg: vsetscscevcevesssenscasvecesusescecssasdesveas lease sapecdcasnetsncionsanse sbeitaoecsenns vereesneeaens 5 


Robust body stature; dorsal body tubercles between the limb insertions in cross-section are 
triangular to elliptical and sharply keeled anteriorly; 34-42 paravertebral tubercles; ventral scales 
are juxtaposed and sharply raised ......... ee eseeeeeseeeeeeeeeees G. toyamai Grismer, Ota & Tanaka, 1994 
Slender body stature; dorsal body tubercles between the limb insertions in cross-section are smooth 
or very weakly keeled; 22—29 paravertebral tubercles; ventrals are flat, wide, and subimbricate to 


TMDTICATE: se cecitaieceetis teditts sedis teat vaseaseeetcrree lander teenth dec eas G. splendens (Nakamura & Uéno, 1959) 
Dorsal banding absent or incomplete ..........cccccesseesseeseeeteeeteeseeeeteeesees G. kuroiwae (Namiye, 1912) 
Prominent dorsal pattern consisting of transverse bands between the nape of the neck and the caudal 
CONSHICUION: 256 5caMaiedeeinie tah dear lieA area iesr iia ied Ableed orien ate tee eset nes 6 


17. 


18. 
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No such reddish or pinkish tint in dorsal pattern and iris «0.0.0.0... G. orientalis (Maki, 1931) 
Reddish or pinkish tint to some extent in dorsal pattern and iris 20... eee eceeeeeececeeeeeeeeeeceeeeeeeeeeeenees 


Precloacal pores in males less than 16 (G. yingdeensis QTOUP) ....cccecceecceesseeeesteeeecneeeeeeeeeeessesteeesaes 8 
Precloacal pores: in:males: more than. 16) ...sc.ssscsegsesseesseseacessdeacessacaats daa cots a satid sans edd odanccsScesaessiaaanes 11 


Scales around midbody 121—128; longitudinal dorsal tubercle rows at midbody 16-17 .............06 
soca esac auabaeatees eased vevetatavistsaasusicastsvedsateetes G. gollum Qi, Wang, Grismer, Chen, Lyu & Wang, 2020 
Scales around midbody 99-115; longitudinal dorsal tubercle rows at midbody 20-28 ..........e 9 


Nuchal loop and body bands with small dark blotches; enlarged row of supraorbital tubercles absent; 
trunk of body usually with a longitudinal light vertebral stripe ........ec ce eceeseesececeeeeeeeeeceeeeaeeseeeenees 
PONS italia nach a8aiaes visas a seteaten eines es seuaceae G. varius Qi, Grismer, Lyu, Zhang, Li & Wang, 2020 
Nuchal loop and body bands without dark blotches; enlarged row of supraorbital tubercles present; 
trunk of body without a longitudinal light vertebral stripe 20.0.0... .ceeeeseeseeeceeeeeeeeeceeceeeeseeeeeneeeaeens 10 


. Tubercles between orbits present; gular scales bordering the postmentals 2-4; preorbital scales 5—7 


a Wehahied saat ivane dee Sadeeaed tuna nec sence tetag edheemuata ds Stesveeedansteadaane cat tezevens G. yingdeensis Wang, Yang & Cui, 2010 
Tubercles between orbits absent; gular scales bordering the postmentals 4—6; preorbital scales 7—9 
daveek Ge Cecetedeaseslaad sotaaueteesaid eon Geottet Gamaseaatveiadetneeneeanete G. zhelongi Wang, Jin, Li & Grismer, 2014 


. Body and limbs robust (G. lichtenfelderi SrOUP) .......ccccccecssessscessceseceseeeseeeseecseceseceseceeeeeseesseecseeeseees 12 
Body and limbs splayed gracile (G. /uii SrOUP) .......:esceeceeseeeeeeeceseeseeseeeeeeaecaeeaeeeeeeaecaeceaeeaeeeeeeeeeaeeas 16 
. Precloacal pores in males 37-46 ....... G. bawanglingensis Grismer, Haitao, Orlov & Anajeva, 2002 
Precloacal pores in males fewer than 37 ........cccccccesccesceeseeesceeseceseceseceseeeseeeseecseceaeenseeeseeeseesseeeaeenseees 13 


. Nuchal loop protracted posteriorly (in V-shape); number of body bands 4; axillary pockets deep; 


precloacal pores in males more than 33 0.0... eee G. zhoui Zhou, Wang, Chen & Liang, 2018 
Nuchal loop rounded posteriorly (in U-shape); number of body bands 3; axillary pockets shallow; 
precloacal pores in males fewer than 33 .0.....eeceeeeseesceseseeeeeeeecesecaeeseeeeeeaecacesaeeeeeaeceeceaeeaeeeeeeeeeaeeas 14 


. Eyelid fringe scales 55—70; enlarged row of supraorbital tubercles absent; scales of upper eyelid to 


top of head equal of those on the top of the head in size «0.0.0... G. hainanensis Barbour, 1908 
Eyelid fringe scales 47-58; enlarged row of supraorbital tubercles present; scales of upper eyelid to 
top of head % of those on the top of the head in SiZe oo... ec eeccceseceeeescesseecseceseceaeeeseeeseesaeeeseenaeens 15 


. Lateral spotting on belly absent; preorbital scales 18—21; scales around midbody 109-118; subdigital 


lamellae under the fourth finger 18-19 oo... ee eeecseesseesteesteeteeeteees G. kwanghua Zhu & He, 2020 
Lateral spotting on belly present; preorbital scales 12—18; scales around midbody 117-130; subdigital 
lamellae under the fourth finger 13-17 oo... eeeecceeseesteesteeeteees G. lichtenfelderi (Mocquard, 1897) 
. Internasal absent .........eeeeeeeees G. catbaensis Ziegler, Nguyen, Schmitz, Stenke & Résler, 2008 
Tnrtemmiasal: present cass cvesuceexsovecsadavets cnssied canes teasuv becca caun chased cxasuegensvpedenass cad eesoscetesaFengesaviedsn dos cesnaeoeaes 17 
Postcloacal tubercles 3-6 wo... cceeecccscceseesessseseeeceseessaeees G. araneus Grismer, Viets & Boyle, 1999 
Postelodcall tubercles 13 cs ssciisassvecatavs dc sddesvecasadseed sageaes casneeacadaneestaasbeenlesaecasadoacdevdenlccuaeeedtteaseeeasocese 18 


Adult body length (SVL) > 126 mm; dorsum without dark blotches .........ccccccsceeseeeteceteeeeeeeeeeeeees 
gaeuitevagedtavisedsinemeeesdmaamtbam srescuesterueten ton G. huuliensis Orlov, Ryabov, Nguyen, Nguyen & Ho, 2008 
Maximum body length (SVL) < 126 mm; dorsum with dark blotches 2.0... ce eceeeeeeseeseeeeeereeeens 19 
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19. Dorsum with many small dark blotches ........ ec eeeeeeeeeeeeees G. luii Grismer, Viets & Boyle, 1999 
— Dorsum with rarely scattered dark blotches... ccceccesceseeeeeeeeeseceaeceseceeceeeeeseesseeaecnseenteseeeeaes 20 


20. Outer surface of upper eyelid composed of granular scales, about the same size of those on top of 


head and with enlarged tubercles ..........cccccecscsseesecseceeeceeceteeeeeeseeeaees G. gezhi Zhu, He & Li, 2020 
— Enlarged row of supraorbital tubercles present; scales of upper eyelid to top of head 4 of those on 
the Top OF the WEAad SIZES: ss sciesis,ckcsess tides veaneeeeecaceged seddedsotandes Seuseus tasted bonacyeesansavens suad@aanaevatbacaevansteaes 21 
21. Precloacal pores in males 31-33 oo... cceccecseesesseesseeseeeteeeteeeees G. kwangsiensis Yang & Chan, 2015 
— Precloacal pores in males fewer than 31... ..ceecesesccsssesceseceeceeeeeeeeseeceaecaceeeseneeaecaeeaeeeneseeceenaeenees 22 


22. Lateral spotting on belly present; nasal scales 6—7; paravertebral tubercles between limb insertions 


30-34; precloacal pores 26-28 oo... ecceeccessceseeteeesseeteeeteeeteeees G. kadoorieorum Yang & Chan, 2015 
— Lateral spotting on belly absent, nasal scales 8-9; paravertebral tubercles between limb insertions 

27-28; precloacal pores 23 o.....ecceceeceeseeeeceeeseeeeeeeeeeaeens G. liboensis Wang, Yang & Grismer, 2013 
Discussion 


Taxonomic review 


Our phylogenetic analyses strongly suggest that populations recorded from Vietnam indeed belong to 
four known species, namely G. catbaensis, G. huuliensis, G. lichtenfelderi and G. luii (Grismer et al. 
1999; Vu et al. 2006; Orlov et al. 2008; Ziegler et al. 2008; Nguyen et al. 2009; Nguyen 2011). Despite 
extensive surveys over the last two decades, G. araneus has not been recorded in Vietnam, leaving 
its accurate type locality in Cao Bang Province, northern Vietnam ambiguous (Grismer et al. 1999; 
Ngo et al. 2016). It is possible that the species has been extirpated from Vietnam as a result of over- 
exploitation to supply the international pet trade or that the species has never occurred at the documented 
type locality (Grismer et al. 1999; Ngo et al. 2016). Similarly, G. /uii has not been recorded again at its 
type locality in China (Grismer et al. 1999; Stuart et al. 2006; Yang & Chan 2015). 


Molecular results supported by our study confirm that G. murphyi, described by Orlov & Darevsky (1999), 
is a junior synonym of G. lichtenfelderi (Grismer 2000). All mainland populations of G. lichtenfelderi 
found in granitic forests are conspecific with those from the granitic offshore islands in Bai Tu Long 
National Park, Quang Ninh Province, northern Vietnam (type locality), with an intraspecific genetic 
distance of less than 0.7% (Supp. file 1: Table S1). However, a few issues need to be resolved in the future. 
In particular, G. hainanensis was recovered as polyphyletic by Liang et a/. (2018) and our phylogenetic 
analyses. Moreover, Zhu et al. (2020a) suggested that G. /uii and G. kadoorieorum are polyphyletic. We 
highly recommend that they may be synonymized. Samples assigned to G. kuroiwae were not recovered 
as monophyletic in our phylogenetic analysis. It is possible that more cryptic species from the Japanese 
group will be discovered. Further studies with additional samples, especially from type localities, should 
be undertaken to clarify these taxonomic problems. 


Regarding morphological analyses, six morphometric characters, namely snout-vent length, head 
length, head width, mouth width, body height and body width strongly supported the overall difference 
among the four species in Vietnam and five characteristics, namely fourth toe, body and limb shapes, 
body bands, precloacal pores and granular scales surrounding dorsal tubercles, mainly accounted for the 
variation among the 17 tiger geckos from China and Vietnam. Six species of the G. kuroiwae group from 
Japan were not included in this study. However, the Japanese G. kuroiwae group can be distinguished 
from Chinese and Vietnamese species by the absence of precloacal pores and unsheathed claws (Wang 
et al. 2014; Yang & Chan 2015; Honda & Ota 2017). We highly recommend that all these traits be 
considered as diagnostic characters for species of Goniurosaurus. 
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Fig. 7. A few key characters for species identification in the genus Goniurosaurus Grismer, Viets & 
Boyle, 1999. A. Snout tip of G. catbaensis Ziegler, Nguyen, Schmitz, Stenke & Résler, 2008 with 
lacking postrostral (internasal) and two supranasals are in contact with each other. B. Snout tip of 
G. lichtenfelderi (Mocquard, 1897) with few internasals. C. Nuchal loop of G. huuliensis Orlov, Ryabov, 
Nguyen, Nguyen & Ho, 2008 protracted posteriorly (in V-shape). D. Nuchal loop of G. lichtenfelderi 
rounded posteriorly (in U-shape). E. Precloacal pores and large swollen hemipenial bulges of G. /wii 
Grismer, Viets & Boyle, 1999 in males, and two postcloacal tubercles. F. Indistinct precloacal pores 
of G. lichtenfelderi in females and only one postcloacal tubercle on each site. Yellow arrows point to 
postcloacal tubercles, and red arrows to other scales. 
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Based on our identification key in accordance with the morphological and phylogenetic analyses, all tiger 
gecko species were grouped into four separate species groups, in agreement with the previous findings 
(Liang et al. 2018; Qi et al. 2020a, 2020b; Zhu et al. 2020a, 2020b). Liang et al. (2018), Qi et al. (2020a, 
2020b) and our genetic analyses revealed that G. bawanglingensis and G. zhoui are embedded within 
the G. lichtenfelderi group, while they were previously placed as sister taxa in the G. /uii group due to 
superficial morphological similarities as a result of adaptation to karst habitat (Grismer ef a/. 2002; Zhou 
et al. 2018). In accordance with these genetic findings, we found that the two karst-adapted species are 
morphologically most similar to three remaining granite-stream-adapted species of the G. lichtenfelderi 
group. Further comparative studies in southern China and Japan may help to resolve the phylogenetic 
relations of the genus Goniurosaurus. 


Due to the high degrees of adaptation to specific microhabitats and local endemism, the genus 
Goniurosaurus could serve as a model system to study the evolution in lizards. In Vietnam, all species 
of Goniurosaurus have been found in isolated geographic ranges and none of them occurs in sympatry. 
For example, the distribution of G. /uii stretches to the north of the Lang Son and Cao Bang provinces 
(Ngo et al. 2016), while G. huuliensis has been recorded in karst forests in the south of Lang Son 
Province. Goniurosaurus catbaensis is currently known only from offshore islands in Cat Ba National 
Park and the Ha Long Bay. Goniurosaurus lichtenfelderi can only be found in granitic habitats at 
several mainland localities and on some offshore islands in the Bai Tu Long Archipelago, which is 
contiguous with Ha Long Bay in the Gulf of Tonkin (Orlov et al. 2008; Gawor et al. 2016; Ngo et al. 
2019a). Goniurosaurus araneus has not yet been re-discovered at its type locality in Cao Bang Province, 
and we only recorded G. /uii inhabiting limestone karst forests of Cao Bang Province (Grismer et al. 
1999; Vu et al. 2006; Ngo et al. 2016). The allopatric distribution has also been documented among 
species of Goniurosaurus from China. To date, neither the four species from the Hainan Archipelago, 
nor the five species of the G. /uii group from Guangxi Province, China, nor the four known species 
of the G. yingdeensis group from Guangdong Province were recorded to occur in sympatry (Grismer 
et al. 2002; Zhou et al. 2018; Zhu et al. 2020a, 2020b). A large river (Zuojiang River) is regarded 
as the potential geographic barrier between the two species G. Juii and G. araneus within Nonggang 
Nature Reserve, China (Chen et al. 2014). Likewise, geographic barriers (e.g., rivers and canyons) 
are considered to prevent a genetic exchange between members of the G. yingdeensis group (Qi et al. 
2020a). The global rising of the sea level during the last melting period of glacial ice has shaped the 
myriad of archipelagoes (Clements et al. 2006; Ziegler et al. 2008; Liang et al. 2018; Ngo et al. 2019a). 
Thus, the oceanic barrier known as an important geographic feature has constrained current distributions 
and limits the genetic exchange of insular tiger gecko species, such as G. catbaensis, species from 
the Hainan Archipelago and G. kuroiwae group (Ziegler et al. 2008; Honda & Ota 2017; Liang et al. 
2018; Zhu et al. 2020b). However, the last glacial maximum period occurring around 10000-50000 
years ago was not sufficiently long for speciation events within Goniurosaurus (Clements et al. 2006; 
Sterling et al. 2006; Ziegler et al. 2008). Indeed, by employing a molecular dating method, Liang et al. 
(2018) estimated that the latest speciation events among tiger geckos took place approximately 2.8—2.9 
Ma, in the late Pliocene, between G. yamashinae-G. orientalis, and G. huuliensis—G. luii. It is likely 
that the global warming during the Pliocene promoted the diversification in the genus, as also reported 
in the crocodile newts (genus Tylototriton Anderson, 1871) in Asia (Bernardes et al. 2020). Thus, the 
genus Goniurosaurus might offer a unique opportunity to investigate mechanisms of speciation and 
evolutionary adaptation, that are expected to highly correlate with a pre-dominantly allopatric mode of 
diversification. 


Implications for conservation 


The taxonomy of several existing species complexes requires further investigation and new species 
are still being discovered. Nevertheless, wild populations of the genus Goniurosaurus are under high 
risks of extinction, due to habitat destruction and exploitation for national and international trade. Of 
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the 23 described tiger gecko species, 12 have been considered threatened (listed as VU, EN and CR) 
and one data deficient (DD) in the IUCN Red List (Ngo et al. 2019b). Assessments for the remaining 
species are still lacking. In order to better regulate the international trade in the species, all species of 
Goniurosaurus from China and Vietnam were recently listed in CITES Appendix II and the Vietnamese 
species were also protected under the Governmental Decree 06/2019/ND-CP (Group IIB). Trade data 
in the US showed that about half of the number of traded tiger geckos were only assigned to the genus 
level (Ngo et al. 2019b). This fact could be a consequence of the highly similar morphology among 
tiger gecko species and lacking guidance to distinguish the species (Ngo et al. 2019b). Besides, the 
conservation status and national protection status are not assessed equally among species. Therefore, 
the present identification guide will aid to prevent tiger geckos from being traded under a wrong name 
to circumvent legislations. The detailed descriptions of morphology provided in this study can assist 
local, national, and international authorities, such as local rangers, CITES authorities, custom officers, in 
enforcing international regulations. The identification key can also be useful for scientists and breeders 
to accurately identify tiger gecko species. However, morphological identification of captive lineages 
should be treated with reservation, as some breeders have produced hybrids (L. Grismer, pers. comm.), 
and thus molecular identification tools should be consulted as well in these cases. 


After the CITES listing, four new species, namely G. gezhi, G. gollum, G. kwanghua and G. varius were 
recently described from China, which are automatically listed in Appendix II of the Convention (Qi 
et al. 2020a, 2020b; Zhu et al. 2020a, 2020b). As the Japanese species appear to be similarly threatened 
and to prevent a shift towards international trade, Janssen & Shepherd (2019) and Ngo et al. (2019b) 
highly recommended that the Japanese species be also included in the CITES Appendices. During 
the time of writing, the Japanese Ministry for Environment announced to list all six endemic species 
of Goniurosaurus from Japan in CITES Appendix III to prevent over-harvesting of wild animals for 
trafficking activities. As a result, the listing of all Japanese species in Appendix III has come into force 
since 14 February 2021 (CITES Notification No. 2020/068). 
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